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ABSTRACT 

We assess what information HST observations of stellar Lya lines can provide 
on the heliosheath, the region of the heliosphere between the termination shock 
and heliopause. To search for evidence of heliosheath absorption, we conduct a 
systematic inspection of stellar Lya lines reconstructed after correcting for ISM 
absorption (and heliospheric/astrospheric absorption, if present). Most of the 
stellar lines are well centered on the stellar radial velocity, as expected, but the 
three lines of sight with the most downwind orientations relative to the ISM flow 
{x^ Ori, HD 28205, and HD 28568) have significantly blueshifted Lya lines. Since 
it is in downwind directions where heliosheath absorption should be strongest, 
the blueshifts are almost certainly caused by previously undetected heliosheath 
absorption. We make an initial comparison between the heliosheath absorption 
and the predictions of a pair of heliospheric models. A model with a complex 
multi-component treatment of plasma within the heliosphere predicts less ab- 
sorption than a model with a simple single-fluid treatment, which leads to better 
agreement with the data. Finally, we find that nonplanetary energetic neutral 
atom (ENA) fluxes measured by the ASPERA-3 instrument on board Mars Ex- 
press, which have been interpreted as being from the heliosheath, are probably 
too high to be consistent with the relative lack of heliosheath absorption seen by 
HST. This would argue for a local interplanetary source for these ENAs instead 
of a heliosheath source. 
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INTRODUCTION 



Ultraviolet observations of nearby stars from the Hubble Space Telescope (HST) have 
proven unexpectedly to be very useful for studying the outermost heliosphere. High resolu- 
tion spectra of the H I Lya line at 1215.67 A sometimes show signatures of absorption from 
H I in the heliosphere, in cases where t he interstellar abso rption is weak enough to allow 
the heliospheric signal to be detectable (IWood et al.ll2005bl ). Absorption is also sometimes 
observed from the astrosphere of the observed st ar, which has provi ded the first detections 
and measurements of winds from solar- like stars (jWood et al.ll2005al ). 



Figure 1 shows a picture of t he basic heli ospheric structure, which is characterized by 
three sharp boundaries (see, e.g., IZankI Il999l ). The innermost is the termination shock. 



where the solar wind is decelerated to subsonic speeds. Then there is the heliopause, which 
separates the plasma flows of the fully ionized solar wind and partially ionized ISM. Although 
the charged constituents of the ISM cannot penetrate the heliopause, the neutral atoms in the 
ISM flow can. The outermos t boundary is the bow shock, where the ISM flow is decelerated 
to subsonic speeds (see, e.g.. IZanklll999t ). The Voyager 1 satellite, which is traveling close to 
the upwind direction of the ISM flow seen b y the Sun, has rec ently crossed the termination 
shock at a distance of 94 AU from the Sun (IStone et al.ll2005l ). roughly consistent with the 
model predictions in Figure 1. 

For most lines of sight, the Lya absorption from heliospheric H I will be dominated by the 
so-called "hydrogen wall" region between the hehopause and bow shock (Region 3 in Fig. 1), 
where models predict tha t interstellar H I will be heated, compressed, an d dece l erated after 



cross ing the bow shock (IBaranov et al.l Il99ll : iBaranov &: Malamal Il993l . Il995l : IZank et al. 



19961 ). There are currently a total of 8 detections of heliospheric absorption. All but one 
are in upwind directions. To be more precise, if we define 9 to be the angle between the 
observed line of sight and the upwind direction of the ISM flow vector seen by the Sun, all 



but on e of the lines of sight with detected heliospheric absorption have 6 < 75° ( IWood et al. 



2005bl ). This is consistent with model predictions, which demonstrate that the deceleration 
of H I within the hydrogen wall relative to the ISM flow is strongest in upwind directions. 
This helps to separate the heliospheric absorption from the ISM absorption, making the 
heliospheric signal easier to detect. 

Neutrals within the heliosphere interact with other particles almost entirely via charge 
exchange, densities being too low for thermal particle collisions to be important. The pop- 
ulation of hydrogen wall H I that dominates Lya absorption in upwind directions is formed 
by charge exchange between interstellar H I that has passed unimpeded through the bow 
shock and ISM protons that have been heated, decelerated, and compressed at the bow 
shock. However, interstellar neutrals can also penetrate the heliopause and charge exchange 
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with solar wind protons in the hehosheath region in between the termination shock and 
hehopause (Region 2 in Fig. 1). Downwind hnes of sight that look through the tail of the 
heliosphere will have extended path lengths through the heliosheath. The Lya absorption 
in such directions, if it can be detected, will potentially be dominated by the heliosheath 
rather than the hydrogen wall. 

In this paper, we focus on observations and model predictions concerning Lya absorption 
from heliosheath neutrals. There are several reasons why this is useful. One is simply that 
heliosheath absorption probes a different region of the heliosphere than the hydrogen wall 
absorption that dominates nearly all previous detections of heliospheric and astrospheric 
Lya absorption. Another reason is that hydrodynamic models of the heliosphere that use 
a fully kinetic treatment of the neutrals seem to predict a disturbingly large amount of 
heliosheath absorption in downwi nd directions, whi c h may be in conflict wi th nondetections 
of absorption in those directions (jWood et al.l 12000 : llzmodenov et al.ll2002l ). We will take a 
closer look at these downwind HST data to see if any evidence for heliosheath absorption 
can be found. If not, we will derive upper limits for the amount of absorption that can be 
present. 

Finally, observed heliosheath absorption can potentially be used to constrain the flux 
of energetic neutral atoms (ENAs) that should be observable within our solar system. Lo- 
cal ENA fluxes should be dominated by neutrals created by charge exchange within the 
heliosheath that are then directed through the termination shock towards the Sun. Thus, 
locally observed ENAs basically sample the same particle population as the heliosheath ab- 
sorption. In 2008, the Interstellar Boundary Explorer ( IBEX) will be launched, which is 
entirely devoted to studying ENAs (IMcComas et al.ll2005l ). Constraints on ENA fluxes from 
Lya absorption will be valuable in preparation for this mission. Furthermore, the ASPERA-3 
instrument on board the Mars Express satellite has already reported a detection of helio- 
spheric ENAs ( iGalli et al.ll2006l ). We will make an initial assessment as to whether the ENA 
fluxes observed by ASPERA-3 are consistent with heliosheath absorption detections and/or 
upper limits. 



2. SEARCHING FOR EVIDENCE OF HELIOSHEATH ABSORPTION 

2.1. A Sirius Example of Heliosheath Absorption 

Currently the only detection of heliospheric Lya absorption in a downwind direction, 
where heliosheath absorption should dominate the hydrogen wall absorption, is for the Sirius 
line of sight at ^ = 139°. The Sirius data are shown in Figure 2. The ISM Lya absorption 
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for this line of s ight is derived with the assistance of constraints provided by other ISM 
absorption hnes (IBertin et al.lll995bl ). There is excess absorption on both sides of the Lya 
hne that cannot be accounted for by the ISM. The excess on the blue side of the line (at 
< —10 km s~^ in the figure) is believed to be from the stellar wind of Sirius, which is a hot 
AO V star (IBertin et al.lll995al : iHebrard et al.lll999l ). The excess on the red sid e of the line 
(at > 40 km s~^), based on the ISM parameters derived by IBertin et al.l (Il995bl ). is claimed 
by llzmodenov et al.l (119991 ) to be heliospheric absorption. 



There are serious problems with this Sirius data, and with the detection of heliospheric 
absorption. The spectral resolution is lower than any other spectrum with a claimed he- 
liospheric absorption detection. Thus, the wavelength calibration will be more uncertain. 
The signal-to-noise (S/N) of the data is also low, as Figure 2 clearly shows. The presence of 
the wind ab sorption complicates the analysis of the blue side of the line, including the D I 



absorption (IHebrard et al.lll999l ). Finally, observations of narrower ISM lines demonstrate 



that there are actually two distinct I SM velocity compon ents for this line of sight, which 



further complicates the Lya analysis. IHebrard et al.l (119991 ) are able to fit HST Lya spectra 



of both Sirius and its companion Sirius B without any need for the existence of excess red- 
side absorption, but this interpretation of the data ideally requires the assumption of a very 
low D/H ratio for the weaker of the two ISM components seen towards Sirius, which would 
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Wood et al. 
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D/H ^ 1.5 X 10"^ 

the presence of excess red-side absorption, possibly heliospheric in origin, to be a plausible 
interpretation of the data. Despite its problems, the line of sight does have the advantage of 
having one of the lowest ISM H I column densities ever observed, due mostly to its very short 
distance of d = 2.6 pc. The IBertin et al.l (Il995bl) analys is suggests an ISM column density 



in cm ^) of logA^n = 17.53, while IHebrard et al.l (119991 ) infer logA^H = 17.81. Both values 



are lower than any other HST-observed downwind line of sight. With the ISM absorption 
being therefore weaker than any other downwind line of sight, this would explain why Sirius 
is so far the only such line of sight with a heliospheric absorption detection. 

Figure 2 also shows the absorption predicted by a hydro dynamic model of th e helio - 
sphere. The model used here is of the type described by iBaranov fc Malamal (Il993l . Il995l ). 
which includes a fully kinetic treatment of the neutrals within the heliosphere. The neu- 
tral hydrogen density, proton density, velocity, and temperature of the undisturbed ISM 
assumed in this model (H I) = 0.18 cm ^, noo(H+) = 0.06 cm ^, Voo = 26.4 km s ^, 

and Too = 6400 K. The top panel of Figure 2 shows absorption predicted for both the up- 
wind 36 Oph line of sight (at 6 = 12°) and the downwind Sirius line of sight {6 = 139°). We 
choose the 36 Oph line of sight for comparison because it represents the most upwind de- 



tection of heliospheric absorption, while Sirius is the only downwind detection (IWood et al. 
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2005bl ). Both absorption profiles are significantly redshifted relative t o the local ISM flow 
speed in those directions, ba sed on the Local Interstellar Cloud vector (iLallement fc Bertin 
I992I : ILallement et al.lll995l ). This is why heliospheric absorption is always found on the 
red side of the ISM absorption. Since astrospheric absorption is observed from outside the 
astrosphere rather than inside, it is conversely found on the blue side of the ISM absorption. 

The absorption predicted for the Sirius line of sight is added to the ISM absorption and 
the result is compared with the data in Figure 2. The model appears to predict too much 
absorption at > 60 km s~^. Similar results h ave previously been found when cornparin g other 
kinetic models with downwind observations (IWood et al.ll200Clt llzmodenov et al.ll2002l ). Since 
charge exchange processes drive the neutrals in the heliosphere out of thermal equilibrium, 
models with a fully kinetic treatment of the neutrals would seem to be necessary to provide 
the most precise description of neutral velocity distributions within the heliosphere, which in 
turn should provide the most accurate Lya absorption profiles. Nevertheless, these kinetic 
models generally seem to overpredict heliosheath absorption in downwind directions. 

However, the nature of the heliosheath absorption predicted by the models is such that 
it is not easy to entirely rule out its presence for observed downwind lines of sight. The top 
panel of Figure 2 shows that the heliosheath absorption towards Sirius is characterized by 
a very broad, extended wing on the red side of the profile. This is very different from the 
absorption profile for the upwind 36 Oph line of sight, which is dominated by the hydrogen 
wall. Heliosheath neutrals are generally much hotter than hydrogen wall neutrals, but he- 
liosheath neutral column densities are much lower. This is why the heliosheath absorption 
tends to be broad but shallow compared to the hydrogen wall absorption. Note that there 
is a weak extended wing for the 36 Oph absorption that would correspond to heliosheath 
absorption in that direction, but the absorption is too weak to be detectable. The weakness 
of the wing is simply because the heliosheath is much narrower in upwind directions (see 
Fig. 1). The hydrogen wall absorption produces a rather sharp, well-defined absorption edge, 
which is something that can be compared with the data in a reasonably definitive manner 
to see whether the model is consistent with the data. However, the broad extended wing of 
the heliosheath absorption in downwind directions complicates matters, for reasons we now 
describe. 

One of the primary sources of systematic uncertai nty in any of these Lya analyses is 
the reconstruction of the background stellar Lya profile (IWood et al.ll2005bl ). No reasonable 
change to an assumed stellar profile will greatly change the velocity where the Lya absorption 
predicted by a model becomes saturated. This is why the sharp absorption edge provided by 
the hydrogen wall absorption provides a much better diagnostic to compare with the data 
than the broad absorption wing that characterizes the heliosheath absorption, which results 
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in the apparent disagreement with the data in Figure 2. The red-wing absorption discrepancy 
in Figure 2 could in principle be resolved by increasing by ~ 30% the background Sirius Lya 
fluxes from 60 — 120 km s~^. Other examples of h ow stellar Lya profi les c an be modified to 



allow such broad wing absorption are presented by lWood et al.l (120001 ) and llzmodenov et al. 



(120021 ). There is still a limit, however, to how much one can modify the stellar profile in 
this fashion before the resulting profile becomes implausible, and there is therefore a limit 
to how much extended wing absorption can be present. For example, modifying the Sirius 
Lya background as described above would create a dubious asymmetry in the stellar Lya 
absorption profile. Quantifying exactly when "dubious" becomes "completely unreasonable" 
is not easy, but that is the task that we now undertake. 



2.2. Stellar Lya Line Bisector Analysis 

Most HST-observed nearby stars that are appropriate targets to se arch for heliospheri c 



and/or astrospheric absorption are cool main sequence stars like the Sun (jWood et al.ll2005bl ) . 
The Lya lines of such stars are relatively narrow, isolated emission lines. Note that Sirius is a 
hot AO V star where the stellar profile is instead a very broad absorption line, so Sirius is un- 
fortunately not a representative example at all. As shown in Figure 2, heliosheath absorption 
will contaminate only the red side of a stellar emission line. Therefore, if there is unrecog- 
nized and undetected heliosheath absorption, the stellar Lya emission profile reconstructed 
from the HST data will end up blueshifted with respect to the stellar rest frame. 

The wings of the Lya line will be formed in the chromosphere of the observed star. 
Although emission lines formed at higher temperatures in the transition re gion can be shifted 



from the stellar re st frame, the wings of Lya should be centered on the star (lAyres et al.lll988 



Wood et al.lll997l ). We can search for previously undetected absorption from the heliosheath 
by searching for cases in which reconstructed Lya profiles end up blueshifted with respect 
to the star. The requirement that a stellar emission line be centered on the stellar radial 
velocity also provides a way to limit modifications to the stellar profile, thereby yielding 
upper limits to the amount of broad heliosheath absorption that can be present in the red 
wing of the Lya line. 



Wood et al.l (j2005bl ) provide a complete list of HST Lya spectra of nearby stars, all 
62 of which have been analyzed to measure ISM H I and D I column densities, search for 
heliospheric/astrospheric absorption, measure chromospheric Lya line fiuxes corrected for 
ISM absorption, etc. From this list we choose a selection of stars to look for shifts in the 
reconstructed Lya line relative to the stellar radial velocity. We exclude hot stars, which have 
Lya in absorption rather than emission. We exclude giant stars, which have Lya emission 
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lines that are too broad to be effective for our purposes. We also exclude short-period, 
unresolved binaries, partly because the stellar rest frame can change due to orbital motion 
during the course of the observation, and partly because in some cases both stars contribute 
to the Lya emission, meaning there is no single stellar rest frame. We only consider data 
with reasonably high S/N, and with a spectral resolution of at least A/AA = 40,000 to 
maximize our confidence in the wavelength calibration. 

And finally, we only consider lines of sight with ISM column densities of logiV(H I) < 
18.3. High ISM columns potentially degrade the effectiveness of our analysis in a couple 
ways. Broad ISM absorption limits us to the far wings of the Lya profile in our search for 
Lya line shifts. Also, the damping wings of the ISM Lya absorption are much stronger 
at higher columns and when the absorption is not exactly centered on the rest frame of 
the star, uncertainties in the ISM H I column density can potentially lead to a shift in 
the reconstructed stellar Lya line that could mimic the shifts we are looking for due to 
heliosheath absorption. 

Our selection criteria reduce the original sample of Lya spectra to 28 targets. These 
are listed in Table 1 in order of increasing 9. Our analysis requires accurate knowledge of 
the stellar photospheric radial velocity, Vrad- The fourth column of Table 1 lists Vrad for 
all of our stars, and the last column provides the references for these values. The quoted 
uncertainties are generally from these references, but we set a floor of 0.1 km s~^ for the 
uncertainties in Table 1. In cases where the observed star is a member of a binary system 
(particularly a Cen B and Ori), we have to consider the effects of orbital motion to 
establish the precise radial velocity at the time of observation, so we have increased the 
quoted uncertainty accordingly. 

Figure 3 illustrates how we measure the velocities of the stellar Lya lines, using eight 
of the stars in Table 1 as examples. The figure shows the observed Lya s pectra and the 



reconstructed stellar line profiles. We refer the reader to IWood et al.l (l2005bl ) and references 
therein for details on the reconstruction process, and for complete illustrations of the line 
profiles. We note that all the stellar profiles used here yield D/H values consistent with 
the D/H = (1.56 ± 0.04) x 10^^ value now believed to apply throughout the L ocal Bubble 



(IWood et al.ll2004l ). Many of the analyses actually assumed this value (e.g., IWood et al. 



2005bh . in which case the H I column density is directly inferred from the more easily 



measured D I column. Since the wings of the H I absorption depend solely on the H I 
column density and not on the doppler parameter, the wings of the stellar Lya profile can 
then be inferred simply by computing the Lya opacity profile, r^, from the H I column 
density and multiplying the observed wings by exprA. 



In Figure 3, we focus only on the lower flux part of the Lya profile where we can measure 
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meaningful line centroids for comparison with the radial velocity. The top parts of the profile 
are not useful for our purposes for two reasons. The first is that the section of the profile 
above the saturated core of the ISM H I absorption is in reality completely unconstrained 
by the data. The second is that cool star Lya line profiles (including t hat of the Sun) can 



have asymmetric self- reversals near line center (see IWood et al.l l2005bl ) , meaning that this 



part of the profile will not necessarily be centered on the star. 

The horizontal dotted lines in Figure 3 indicate flux levels where we measure the centroid 
of the line. Collectively these centroids form a line bisector. These bisectors are shown as 
near-vertical solid lines in the figure. For five of the stars (70 Oph A, e Ind, ( Dor, r Get, 
and e Eri) the bisector ends up where we expect it, very close to the stellar radial velocity. 
However, for the three most downwind lines of sight {x^ Ori, HD 28205, and HD 28568) the 
bisector is significantly blueshifted with respect to the star. This is exactly the signature of 
undetected heliosheath absorption that we are looking for, and the downwind direction is 
where we expect to see it. 

We compute the average velocity of the bisector, Vhis, for each star, ignoring the lowest 
two flux levels where low S/N in the extreme far wings of the line sometimes causes unrea- 
sonable excursions in the bisector (e.g., HD 28568 in Fig. 3). These bisector velocities are 
listed in Table 1. We define the discrepancy between these velocities and the stellar rest 
frame as AV = Vus — Vrad- In Figure 4, AV is plotted versus 9 for all the stars listed in 
Table 1. 

For the three 6 > 160° cases the Lya profile is significantly blueshifted. This blueshift 
is common to all 6 > 160° lines of sight, is not seen in other directions, and is exactly in the 
direction where we would expect heliosheath absorption to be maximized. Thus, we conclude 
that it is highly likely that heliosheath absorption is responsible for these Lya blueshifts. 
For each of these stars, we can construct a stellar profile that is forced to be centered on the 
stellar rest frame by reflecting the blue wing of the original reconstructed profile across the 
stellar radial velocity onto the red wing. The flux difference between the red wing of this 
profile and the original one provides an estimate of the heliosheath absorption that is present, 
which can be compared with the predictions of hydrodjTiamic models of the heliosphere. 

Other than the three 6 > 160° sight lines, the reconstructed Lja profiles are encour- 
agingly close to their expected locations at the stellar radial velocity. The average and 
standard deviation of the 6 < 160° data points is AV = 0.0 ±1.4 km s~^. The scatter can 
be easily explained by uncertainties in the spectral wavelength calibration, and systematic 
errors present in the profile reconstruction process. We note that most of the HST data for 
the stars listed in Table 1 are STIS/E140M spectr a, which have demo nstrated la random 



wavelength calibration uncertainties of ±0.9 km s ^ ( IWood et al.ll2005bl ). This same analysis 
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also found a systematic offset of —1.2 km s~^ for STIS/E140M Lja spectra. We corrected all 
the STIS/E140M spectra for this shift before starting the bisector analysis. The centering of 
the AV values on km s~^ (for 6 < 160°) implies that this correction successfully removed 
all significant systematic wavelength calibration problems. 

None of the 6 < 160° profiles is shifted away from the stellar rest frame by more than 3 
km s~^. These results imply that a reconstructed stellar Lya line should be within 3 km s~^ 
of the stellar rest frame to be reasonable. This limits how much the Lya profiles can be 
changed to allow for potential heliosheath absorption. We can therefore use these limits to 
establish upper bounds on the amount of heliosheath absorption that can be present for the 
6 < 160° lines of sight listed in Table 1. 

To construct an upper limit to the red wing of the Lya profile, we define an upper limit 
bisector that at all flux levels is +3 km s~^ from the stellar radial velocity or +3 km s~^ from 
the velocity of the original bisector, whichever is largest. An upper limit for the red wing 
of the Lya profile is then derived by reflecting the blue wing of the original reconstructed 
profile across the upper limit bisector. Figure 5 shows the result for the case of r Get. The 
difference between the upper limit profile and the original one represents an estimate of the 
maximum amount of heliosheath absorption that can be present in that direction. 

One final issue worthy of mention concerns astrospheric absorption. If we observed a 
star through the downwind tail of its astrosphere, it is possible that astrosheath absorption 
would lead to a reconstructed stellar Lya profile that is redshifted relative to the star rather 
than blueshifted, keeping in mind that astrospheric absorption appears on the opposite side 
of the Lya line from the heliospheric absorption (see §2.1). However, Figure 4 shows no data 
points with AV > +3 km s~^. Given that a heliospheric effect is only seen for 6 > 160°, we 
conclude that either none of these 28 sight lines has 6 > 160° for the stellar astrosphere, or 
if there are any 6 > 160° astrospheric lines of sight the ISM surrounding the star contains 
no neutrals to provide Lya a bsorption, which sh ould be a common occurrence within the 
mostly ionized Local Bubble (IWood et al.ll2005bl ). 



3. MODELS WITH SINGLE- AND MULTI-COMPONENT PLASMAS 



As an example of how the heliosheath absorption detections and upper limits can be 
used, we now compare the Lya data with the predictions of a pair of heliospheric models. 
One of t he models is the kinetic model de scribed in §2.1 and used in Figure 2. This model 
is of the iBaranov fc Malama] (jl993l . Il995l ) type, which uses a fully kinetic treatment of the 
neutrals but treats the plasma as a single hydrodynamic fluid. This plasma treatment is 
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simplistic. The best evidence for this concerns pick-up ions inside the termination shock 
(Region 1 in Fig. 1). Pick-up ions originate from charge exchange between outflowing solar 
wind protons and inflowing ISM neutrals, which creates a solar wind neutral H atom and an 
inflowing ion that is then "picked-up" by the solar wind. Observations show that the pick-up 
ions are not thermalized wi th the solar wind, and th e pickup ions also show non-Maxellian 



velocity distributions (e.g.. iGloeckler fc Geisd l2004l ). Thus, combining solar wind plasma 



and pick-up ions into a single fluid in a heliospheric model, as is typically done, will lead to 
inaccuracies in the particle distributions. 



Malama et al.l (120061 ) have recently attempted to solve this problem by modifying the 
Baranov Sz Malama heliospheric modeling code to allow a complex multi-component treat- 
ment of the plasma, as well as including the usual full kinetic treatment of the neutrals. 
We consider here a multi-component model of this type using the same input parameters 
assumed in the single-component model (see §2.1). Figure 1 shows the differences in global 
heliospheric structure suggested by the two models. 

Figure 6 compares the absorption predicted by the single- and multi-component models 
with the Lya absorption observed for four downwind lines of sight, including the Sirius line 
of sight discussed in §2.1. The figure focuses on the red side of the absorption where the 
heliospheric absorption should be located. Using the original reconstructed stellar profiles, 
the ISM absorption alone is illustrated, and the absorption predicted by the heliospheric 
models is shown after being added to the ISM absorption. 

Dashed lines in Figure 6 indicate the extra heliosheath absorption that is allowable 
by the data based on the upper limit stellar profiles derived as described in the previous 
section (see Fig. 5). Models can be considered to be consistent with the data if they do not 
predict more absorption than these dashed lines, so the dot-dashed lines should lie above the 
dashed lines in Figure 6. The dashed lines have a lower velocity limit in Figure 6. This is 
because the line bisector analysis described above is not applicable close to line center (for 
reasons discussed in §2.2), so we cannot really quantify an upper limit to the profile near 
line center. No absorption limits have been derived for Sirius since it is a hot star where 
the stellar profile is a broad absorption line rather than a comparatively narrow emission 
line (see §2.2). However, as described in §2.1, this is the one downwind line of sight with a 
previous detection of heliosheath absorption, with the ISM absorption underpredicting the 
observed absorption in Figures 2 and 6. 

As described in §2.2, the Lya bisector analysis suggests that O^i and two other 
very downwind lines of sight (HD 28205, HD 28568) show evidence for broad heliosheath 
absorption. The thick solid line in the panel of Figure 6 corresponds to the amount 

of excess absorption that should be present to allow a stellar line profile that is centered on 
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the rest frame of the star. For Oii, therefore, the model absorption in Figure 6 is best 
compared with the thick sohd hne rather than with the data. We could have added panels 
for HD 28205 and HD 28568 to Figure 6, but this would be somewhat redundant since the 
Ori, HD 28205, and HD 28568 lines of sight are very similar, and the x^ O^'i data are in 
any case far superior to the others in terms of S/N. The negative radial velocity of x^ Ori 
profile also makes that line of sight superior to HD 28205 and HD 28568. A blueshifted 
profile means that the bisector method can detect and quantify heliosheath absorption at 
lower heliocentric velocities than for redshifted profiles like those of HD 28205 and HD 28568, 
since you cannot do the search close to the center of the stellar profile. For x^ On, we can 
extend the dashed and solid lines down to ~ 80 km s~^ (see Fig. 6), but for HD 28205 and 
HD 28568 this could only be done down to ~ 130 km s~^. 

The absorption limits (dashed lines in Fig. 6) allow a significant amount of broad he- 
liosheath absorption to exist while still being consistent with the data. Thus, the initial 
impression of disagreement between the models and the data in Figures 2 and 6 is partly 
illusory. The multi-component model predicts less absorption than the single-component 
model. This provides better agreement with the data, particularly for the Sirius line of 
sight. The absorption predicted by the single-component model may be in conflict with 
the absorption limit towards e Eri near 100 km s~^, but the conflict is greatly lessened for 
the multi-component model. The greater amount of broad heliosheath absorption for the 
single-component model indicates that the simplistic treatment of the plasma in this model 
overemphasizes the effects of minority high-temperature constituents on the plasma velocity 
distribution. The effect is to artificially broaden the distribution, a problem that is then 
transmitted to the neutral distributions via charge exchange. 

One very important caveat that must be mentioned is that the heliospheric models that 
we have employed here use a spatial grid that extends only 700 AU from the Sun (see Fig. 1), 
which may not extend far enough to capture all the heliospheric absorption in downwind 
directions along the lengthy tail of the heliosphere. Thus, the model absorption predictions 
in Figure 6 may end up being lower limits to the actual amount of absorption that these 
models would really predict given a sufficiently large grid. This problem could potentially 
worsen agreement between the models and the data, particularly for the most downwind 
directions where the bisector analysis has provided new evidence for heliosheath absorption. 
In the future, we will expand the grid for the mult i- component m odel up to 10,000 AU, as 



has already been done for sin gle-fluid Baranov & Malama models ( llzmodenov fc Alexashov 



200a lAleksashov et al.ll2004f ). 
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CONNECTING HELIOSHEATH ABSORPTION WITH LOCAL ENA 

FLUXES 



The ASPERA-3 instrument on Mar s Express was designed prim arily to study particles 
escaping from the Martian atmosphere (IBarabash fc Lundinl |2006| ). However, ASPERA-3 
has detected ENAs even when not po inted at Mars, both en route to Mars and while in orbit 
around the planet. iGalli et al.l (120061 ) and lWurz et al.l (120061 ) have argued that a heliosheath 
origin for these particles is more likely than a more local interplanetary source. 

The particle energy spectra observed for th e nonplanetary E NAs are nearly always 
consistent with the following power law behavior ( iGalli et al.ll2006l ): 



^-1.6 
^-3.3 



for 
for 



E < 0.77 keV 
E > 0.77 keV. 



:i^ 



Particle fluxes are variable, though the nature of the ASPERA-3 observations makes it 
difficult to distinguish temporal from spatial variability. Within ASPERA-3's energy range 
of 0.2 — 10 keV, logarithmic particle fluxes are in the range log Fas = 3.7 — 5.0 (in cm~^ 
s~^ ster~^ units). If the nonplanetary ENAs observed by ASPERA-3 are indeed heliosheath 
neutrals, they then represent the same particle population that we have been investigating 
with regards to the Lya absorption produced by these neutrals. Thus, we can in principle 
test the heliosheath interpretation of the nonplanetary ASPERA-3 ENAs by seeing if the 
observed fluxes are consistent with the Lya data. 

Dotted lines in Figure 7 illustrate the observed energy spectrum given by equation (1), 
with the lower and upper lines indicating the lower and upper range of observed fluxes. The 
lower bound of ASPERA-3's energy range, 0.2 keV, represents a particle velocity of 196 
km s~^. This corresponds with the upper bound of the 50 — 200 km s~^ velocity range where 
the Lya data provide the best absorption constraints (see Fig. 6). Thus, before we can even 
hope to use the ASPERA-3 measurements to predict Lya absorption, we must extrapolate 
the ASPERA-3 energy spectra to lower energies. There are three different extrapolations 
assumed in Figure 7. The solid line simply continues the E~^'^ power law to lower energies. 
The dashed line assumes a flat spectrum below 0.2 keV, and the dot-dashed fine assumes 
a drop in flux below 0.2 keV. The four Figure 7 panels show the resulting spectra for four 
different values of log Fas- 

The biggest difficulty in connecting the ASPERA-3 ENA and HST Lya data lies in 
the issue of how to extrapolate a locally observed particle spectrum along an entire line 
of sight through the heliosphere. If the heliosheath is the source of the ENAs, then for 
our purposes the region inside the termination shock (Region 1 in Fig. 1) is in effect an 
empty cavity that is being filled by a flux of heliosheath particles from all directions. The 
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heliosphere is not spherically symmetric, so it is likely that fluxes will be larger from some 
directions than others. Nevertheless, assuming that within the termination shock the particle 
distributions are isotropic and the same throughout Region 1 may be a reasonable first-order 
approximation. This will probably not be the case within the heliosheath itself. Thus, for 
our purposes here, we will confine our attention only to absorption by heliosheath neutrals 
within the termination shock. 

In downwind directions, which are the most appropriate for heliosheath Lya absorption, 
the termination shock distance is about 200 AU (see Fig. 1). The inset figures in Figure 7 
show the Lya absorption produced by heliosheath neutrals within this distance, assuming 
that the ENA spectra shown in the figure are isotropic and apply along the entire 200 AU 
path length. The absorption limits derived for the r Get and e Eri lines of sight in §2.2 (see 
Fig. 6) are also shown, for comparison with the absorption predicted by the ENA spectra. 

The two highest flux values explored in Figure 7, logF^s = 4.55 and logF^s = 4.95, 
clearly lead to too much Lya absorption to be consistent with the limits from the HST data. 
A flux of log Faz = 4.15 is only consistent with the Lya data if fluxes decrease below 0.2 keV, 
and even the lowest flux of logF^a = 3.75 is only consistent with the absorption limits if the 
particle fluxes are flat or decrease below 0.2 keV. The situation is actually even worse than 
this, given that the absorption predictions in Figure 7 only include absorption from within 
the termination shock. Considering that a downwind line of sight will have a path length 
through the heliosheath much longer than the ~ 200 AU path length inside the termination 
shock, it is clear that including absorption from the heliosheath could easily increase the 
amount of absorption many times above the predictions shown in Figure 7. With this in 
mind, it is very questionable whether any of the energy spectra in Figure 7 are consistent 
with the relative lack of heliospheric absorption observed in downwind directions. At the 
very least, the spectrum must lie near the lower bound of the ASPERA-3 range, with fluxes 



that d ecrease to lower energies. It is worth noting that the flux range quoted by iGalli et al. 



(120061 ) do es not take into a ccount numerous upper limits when no EN As were detected (see 
Fig. 4 in iGalli et al.l l2006l ). suggesting that perhaps the lower bound could be decreased 
further without being truly discrepant with the ASPERA-3 data. 

Nevertheless, our difficulty in reconciling the ASPERA-3 ENA fluxes with the HST 
Lyg spectr a is a s trong argument against the ASPERA-3 ENAs being heliosheath neutrals. 
Galli et al.l (120061 ) acknowledge other difficulties with the heliosheath interpretation. One is 
that many heliospheric models predict ENA fluxes at Earth about an order of magnitude 
lower than those observed by ASPERA-3 (IGruntman et al.ll200ll : IWurz et al.ll2006l ). The sub- 
stantial variability seen by ASPERA-3 is also potentially difficult for a distant heliosheath 
source to explain, and is more suggestive of a local interplanetary source. Identifying the ori- 
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gin of the ASPERA-3 ENAs is particularly important considering the upcoming 2008 launch 
of IBEX, which unlike ASPERA-3 is entirely designed and devoted to studying heliosheath 
ENAs. If there is an unexpectedly strong local interplanetary ENA source that is responsible 
for the ASPERA-3 ENAs, IBEX will presumably see it as well and will somehow have to 
distinguish between these particles and the true heliosheath ENAs to accomplish its mission. 
It is also worth noting that the ASPERA-4 instrument on board Venus Express will soon be 
able to confirm or refute ASPERA-3's measurements of nonplanetary ENAs. 

Based on Figure 7, we quote a very conservative upper limit of F < 5000 cm^^ s^^ ster^^ 
for the local flux of heliosheath ENAs in the 0.01—0.2 kcV energy range most applicable to the 
Lya absorption data. The ASPERA-3 and ASPERA-4 instruments are not sensitive to these 
low energies, but IBEX will be. Unless the heliosheath unexpectedly produces significantly 
less absorption than the region inside the termination shock, the actual 0.01 — 0.2 keV flux 
can be expected to be well below this limit. More precise constraints on ENA fluxes from 
the Lya data will require the assistance of kinetic heliospheric models to better connect 
heliosheath absorption along a downwind line of sight with the ENA energy spectrum close 
to the Sun. Once IBEX is launched it will be very interesting to see if heliospheric models 
are capable of simultaneously reproducing both the local ENA spectrum and the heliosheath 
absorption data, particularly the previously detected absorption towards Sirius and the three 
new 9 > 160° detections reported here for the first time. 



5. SUMMARY 

The absorption within the Lja lines of nearby stars sometimes has a detectable he- 
liospheric contribution, in addition to the ubiquitous interstellar H I and D I absorption. 
Most detections are in upwind directions, where the hydrogen wall dominates the absorp- 
tion. However, in this paper we have focused on absorption from the heliosheath, which is 
only prominent in downwind directions. Our findings are as follows: 

1. In order to search for evidence of previously undetected heliosheath absorption, we have 
measured the bisectors of 28 stellar Lya lines reconstructed from the data after cor- 
recting for the broad ISM absorption (and heliospheric/astrospheric absorption, when 
present). We find that the stellar Lja profiles of the three most downwind lines of sight 
(X^ Ori, HD 28205, and HD 28568), all with 6 > 160°, are significantly blueshifted rel- 
ative to the stellar radial velocities. This is indicative of broad absorption in the red 
wing of the Lya line that was not accounted for in the reconstruction of the stellar 
profile. This is exactly the signature of heliosheath absorption we are looking for, and 
the most downwind lines of sight are where we expect to see it. Thus, we add these 
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three lines of sight to the hst of eight previous hehospheric Lya absorption detections, 
only one of which (Sirius) is a downwind line of sight. For the three new detections, 
the amount of heliospheric absorption can be estimated by reanalyzing the Lya data 
using a stellar profile that is forced to be centered on the stellar rest frame. 

2. The stellar Lya profiles of the 25 lines of sight with 9 < 160° are all found to be within 

±3 km s^^ of the stellar radial velocity. This represents a conservative quantitative 
estimate of the accuracy of the reconstructed Lya profiles. We can define an upper 
limit to the red wing of the stellar Lya profile by requiring the bisector to be no more 
than +3 km s~^ from either the stellar radial velocity or the bisector of the original 
reconstructed profile, whichever is greatest. In this way, we derive for the first time 
quantitative upper limits for the amount of broad heliosheath absorption that can be 
present in the red wings of Lya lines without clearly detectable heliosheath absorption. 

3. We compare the Lya absorption predicted by two kinetic heliospheric models with down- 

wind Lya lines observe d by HST. One model uses the complex mult i- component plasma 



treatment discussed by lMalama et al.l (120061 ) . while the other uses a simple single-fluid 



approach. The multi-component model predicts significantly less heliosheath absorp- 
tion than the single-fluid model, in better agreement with the data. This suggests 
that the single-fluid plasma approach, which is the usual assumption in heliospheric 
models, may not reproduce particle distribution functions accurately enough for some 
applications. 

4. Meaningful comparisons between models and the Lya absorption data in downwind di- 

rections require the use of a model grid that extends a sufficiently long distance down 
the tail of the heliosphere to capture all of the heliosheath absorption. Most hydro- 
dynamic models of the heliosphere, including the ones used here, do not extend the 
grid far enough for these purposes. In the future, we hope to perform more precise 
data/model comparisons using sufficiently large grids. 

5. The ASPERA-3 instrument on board Mars Express has detected a substantial flux of non- 

planetary ENAs, which have been interpreted as being heliosheath neutrals. We believe 
that the observed ENA fluxes likely predict too much heliosheath Lya absorption in 
downwind directions to be consistent with the HST Lya spectra. This represents a 
strong argument against a heliosheath origin for these ENAs, instead suggesting a local 
interplanetary source. 

6. The HST Lya absorption data can provide useful constraints on the heliosheath ENA 

fluxes that IBEX can expect to observe when launched in 2008. We here quote a very 
conservative upper limit of F < 5000 cm~^ s~^ ster~^ for the local flux of heliosheath 
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ENAs in the 0.01 — 0.2 keV energy range most applicable to the Lya absorption data. 
We hope to compute more stringent constraints in the future using kinetic models 
of the heliosphere to connect particle distributions in the heliosheath with those that 
IBEX will observe locally. 

We would like to thank Peter Wurz for providing information on the ASPERA-3 in- 
strument and its measurements. This work was supported by NASA grant NNG05GD69G 
to the University of Colorado. V. I. was also supported by RFBR grant 04-02-16559, the 
"Dynastia" Foundation, and the "Foundation in Support of Russian Science" . 
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Table 1. Lyct Line Velocities 



Star 


6» 

fdeg) 


(km s^i) 


y ^ 

(km s""'') 


Ref. 


36 Oph A 


12.3 


0.2 


0.5 ±0.2 


1 


70 Oph A 


23.8 


-6.8 


-7.9 ±0.5 


2 


HD 165185 


25.4 


13.5 


14.9 ±0.1 


3 


HD 128987 


33.3 


-22.8 


-23.0 ±0.1 


3 


^ Boo A 


46.2 


0.2 


1.3 ±0.1 


4 


a Cen B 


52.4 


-20.7 


-20.8 ±0.5 


5 


61 Vir 


53.1 


-7.5 


-7.9 ±0.1 


4 


X Her 


59.6 


-54.1 


-56.1 ±0.1 


4 


e Ind 


69.8 


-38.9 


-38.9 ±0.3 


6 


61 Cyg A 


78.3 


-65.3 


-65.7 ±0.1 


4 


HD 116956 


84.4 


-12.8 


-12.3 ±0.1 


3 


HD 97334 


96.6 


-2.0 


-3.7 ±0.1 


4 


EV Lac 


98.2 


3.0 


0.4 ±0.1 


4 


C Dor 


107.7 


-1.4 


-0.7 ±0.1 


7 


HD 166 


111.0 


-7.5 


-6.5 ±0.1 


4 


DK UMa 


111.5 


-27.6 


-27.3 ±0.1 


8 


/3 Cas 


111.8 


13.3 


11.3 ±0.5 


9 


HD 73350 


122.1 


33.8 


35.4 ±0.1 


4 


T Get 


122.9 


-18.2 


-16.6 ±0.1 


4 


HD 43162 


137.3 


22.1 


21.7±0.1 


3 


HD 37394 


140.3 


1.1 


1.2 ±0.1 


4 


€ Eri 


147.4 


17.9 


16.3 ±0.1 


4 


5 Eri 


148.8 


-8.8 


-6.3 ±0.1 


4 


K, Get 


152.5 


16.4 


19.0 ±0.1 


4 


40 Eri A 


154.7 


-42.4 


-42.3 ±0.1 


4 


xi Ori 


166.3 


-21.1 


-15.3 ±0.5 


10,11 


HD 28205 


172.0 


33.2 


39.3 ±0.2 


12 


HD 28568 


173.3 


30.3 


40.9 ± 1.3 


12 



References. — (1) Beavers & Bitter 1986. (2) Tolcovinin 
& Smekhov 2002. (3) Gaidos et al. 2000. (4) Nidever et al. 
2002. (5) Pourbaix et al. 2002. (6) Buscombe & Kennedy 
1968. (7) Murdoch et al. 1993. (8) Duquennoy et al. 1991. 
(9) Duflot et al. 1995. (10) GuUberg & Lindegren 2002. 
(11) Irwin et al. 1992. (12) Ferryman et al. 1998. 
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Fig. 1. — The locations of the termination shock (TS), hehopause (HP), and bow shock 
(BS) according to the single-fluid (solid lines) and multi-component (dashed lines) plasma 
models used in the text (see §2.1 and §3). The Sun is at the origin. The ISM flow direction 
is indicated by the Yism vector, and the angle with respect to the upwind direction of this 
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Velocity (km s ^) 



Fig. 2. — The top panel compares heliospheric H I Lya absorption profiles plotted on a 
heliocentric velocity scale predicted for the upwind 36 Oph (dot-dashed line) and downwind 
Sirius (dashed line) lines of sight, based on a hydrodynamic model of the heliosphere (see 
text). The thick vertical lines indicate the projected ISM flow speed in these two directions. 
The bottom panel shows an HST spectrum of the Lya absorption line observed towards 
Sirius (histogram). The dotted line is the ISM absorption alone for this line of sight. The 
broad absorption centered at +20 km s^^ is interstellar H I absorption and the weak ab- 
sorption at about —65 km s~^ is from interstellar deuterium. The thick solid fine is the 
predicted heliospheric absorption from the upper panel (dashed line) combined with the 
ISM absorption. 
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Velocity (km s ^) 



Fig. 3. — Lya line bisector measurements for various stars. The histograms are the HST 
Lja spectra. The broad absorption dominating the observed hne profiles is H 1 absorption, 
and the weak absorption about —82 km from the center of the H I absorption is from 
interstellar deuterium. The solid lines above the data are the reconstructed stellar Lya 
profiles. The dotted fines indicate the fiux levels where the bisectors of the stellar profiles 
are measured, which are the thick near-vertical solid lines. The vertical dashed fines are 
the stellar radial velocities. Note that for the three most downwind lines of sight (x^ Ori, 
HD 28205, and HD 28568) the bisector is blueshiftcd relative to the stellar rest frame, 
suggesting the presence of heliosheath absorption in the red wing of the line. 
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Fig. 4. — Velocity differences between the Lya line bisector and the stellar rest frame ( = 
His — Kad) plotted as a function of ^, the line-of-sight orientation with respect to the ISM 
flow vector. There is evidence that Lya profiles reconstructed for the most downwind lines 
of sight {Q > 160°) are systematically blueshifted with respect to the stellar rest frame, 
suggesting the presence of heliosheath absorption. 
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Velocity (km s ^) 



Fig. 5. — The reconstructed stellar Lya profile for r Get (solid line), and an upper limit for 
the red wing of the profile (dashed line) derived by requiring that the profile be within 3 
km of the stellar radial velocity (see text). 



-26- 




Fig. 6. — Normalized Lya spectra in downwind directions, focusing on the red side of the 
Lya absorption hne where the hehospheric absorption resides. The dotted hne is the ISM 
absorption alone, which can account for all the observed absorption, except for Sirius. The 
thin and thick dot-dashed lines are the predicted hehospheric Lya absorption of models with 
a single plasma fluid component and a multi-component plasma treatment, respectively (see 
text). The dashed lines (when present) are the maximum heliospheric absorption that is 
allowable by the data, based on upper limit stellar Lya profiles derived as described in the 
text. The models can be considered consistent with the data if their predicted absorption 
lies above these dashed lines. For Ori, there is evidence from a blueshift of the original 
reconstructed Lya line that there is extended heliosheath absorption present, and the thick 
solid line is the amount of absorption that should be observed if the real stellar Lya line is 
centered on the stellar rest frame. 
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Fig. 7. — The dotted lines in each panel encompass the range of ENA spectra observed 
by ASPERA-3, corresponding to different 0.2 — 10 keV particle fluxes (in cm~^ s~^ ster~^ 
units) of log Faz = 3.7 — 5.0. Assuming different fluxes within this range, each panel shows 
three different spectra, which use different extrapolations to low energies. The inset figures 
show the H I Lyct absorption from within the termination shock predicted by these spectra, 
plotted on a heliocentric velocity scale (in km s~^). The two thin solid hnes that stretch 
from ~ 70 to 175 km s~^ are the absorption limits derived for r Get and e Eri based on the 
Lya bisector analysis (sec §2.2 and Fig. 6). The absorption predicted by the energy spectra 
must lie above these limits to be consistent with the Lya data. 



